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Abstract 

We consider here quark matter equation of state in a relativistic harmonic 
confinement model at zero temperature. The same is considered to study 
phase transition of neutron matter to quark matter at high densities. This, 
along with a phenomenological equation of state in the neutron matter sector 
is used to study hybrid stars. Using Tolman Oppenheimer Volkoff equations, 
stable solutions for hybrid stars are obtained with a maximum mass of 1.98 
Mq and radii around 10 kms with a quark core of about 1 to 2 kilometers. 
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Typeset using REVT^ 



I. INTRODUCTION 



It is believed that hadronic matter undergoes a phase transition to quark matter at high 
temperatures and/or high densities. We may have these transitions possible in heavy ion 
collision experiments. However, signals here may get masked since such a phase lasts for a 
very short time followed by hadronisation. On the other hand we could also have compact 
stellar objects formed in a gravitational collapse where quark gluon matter could exist at 
sufficiently high densities. Since quarks and gluons are confined with absence of totally 
asymptotic states the overall confinement can affect the equation of state for quark matter. 
A phase transition of neutron matter to quark matter at zero temperature or temperature 
small compared to degeneracy temperature allows the existence of hybrid stars i.e., stars 
having a quark core and crust of neutron matter with appropriate pressure balancing at the 
interface [1-6]. In fact, the large magnetic field of the pulsars could be held by quark matter 
with nonzero electrically charged constituents rather than neutron matter. 

Although there is a naive wish that quantum chromodynamics (QCD) at high tem- 
perature and/or densities should approach free field limits, there have been many rather 
convincing results that QCD remains as nonperturbative as the confined phase at high tem- 
peratures 0]. Expecting a similar behaviour at finite baryon densities also, we are motivated 
to study thermodynamic properties, e.g., the equation of state at finite densities including 
the effect of confinement explicitly. 

In the present paper we first analyse the quark matter equation of state that explicitly 
takes into account the overall confinement of the quarks using a relativistic harmonic oscil- 
lator potential. The harmonic confinement model as used here has been applied earlier for 
study of hadronic properties as well as to consider equation of state at finite temperature and 
zero baryon density [§. We consider this model at finite baryon density as may be relevant 
for hybrid stars. In the neutron matter sector we shall use a phenomenological equation of 
state. This is then used to consider phase transition to quark matter and to study hybrid 
stars. 
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It is generally felt that it is difficult to take correct thermodynamic limit in the presence 
of confinement. However the oscillator states of the relativistic harmonic confinement model 
are exactly analogous to the Landau orbits of electrons in a magnetic field ||^. The Landau 
orbits are characterised by an arbitrary centre of confinement. The collection of such orbits 
has no difficulty in forming a large system with appropriate thermodynamic limit. However, 
if few such orbits form a colourless bound state like a hadron the centre of mass of the 
system is free to move and the centre of mass correction of the bound state has to be 
applied. In the context of neutron stars the gravitational field crunches the neutrons and 
a large macroscopic quark matter consisting of the Landau like states is formed. In any 
confinement model the energy per particle grows rapidly with the number of particles. This 
would mean that one requires infinite energy in the thermodynamic limit. However, one has 
to keep in mind that the gravitational potential which grows as square of the number of 
particles is primarily responsible in producing the macroscopic state. The thermodynamic 
properties of this system are calculated by using the partition function which does not 
involve gravitational field but takes account of modified phase space occupation because of 
long range confinement. Then the full system is considered inclusive of gravitational field 
in TOV equations. The analogue of an electron in a magnetic field is used for obtaining the 
statistical weights in converting momentum space integrals to sums over possible states in 
the confining field. 

We organise the paper as follows. In section 2, we derive the equation of state for quark 
matter in relativistic harmonic confinement model and use Gibbs criterion to consider the 
possibility of a phase transition of neutron matter to quark matter. In section 3 we shall 
use these equations of state to consider hybrid stars using Tolman Oppenheimer Volkoff 
equations. In section 4, we discuss and summarise our results. 
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II. RELATIVISTIC HARMONIC CONFINEMENT MODEL AND QUARK 

MATTER EQUATION OF STATE 



In the harmonic confinement model the relativistic equation for the quarks is given 

as 

[tYd^-M,-V{r)]^,{r} = (1) 

where V{r) is the 'potential' with a Lorentz scalar plus vector structure given as 

r(r) = l(l + 7o)«V. (2) 

In the above, the subscript 'i' in iplr) represents the colour and flavour index of the quark. 
This model had been applied earlier to study hadronic properties and nucleon nucleon 
angular distributions JlOl- In terms of the two component single particle wave functions, 
equation (1) may be written in the simple form 

[p' + fiV2]0,(f) = (E' - Mf )0,(r-O, (3) 

where Q = a\jE + Mj are "frequencies" (assumed to be constant) of the quark flelds. We 
may note that the confining term fi^r^ is given by Tr{^J^^A^^\ where J^^ is the background 
field similar to quantum motion of an electron in external magnetic field giving rise to 
Landau orbits. The 'potential' V^(r) in equation (2) which has both Lorentz scalar and 
vector components gives rise to the energy levels as H,^ 

= M2 + (2n + l)n (n = l,2,---). (4) 

For the study of hadronic properties the parameters chosen were 

= 160.66Me1/ 

Ms = mMMeV 
Vn = 1.162fm-\ 
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which we shall be using to derive the quark matter equation of state. The thermodynamic 
pressure Pi (i=u,d,s) is then given by [Tl[ 



^ = 9V9I E 9n [log(l + e-f^(^^->^A + log(l + e-^(^"+-'))] , (5) 

where gy, gi and gn are respectively the statistical weights corresponding to volume, internal 
degrees of freedom (spin and colour) and the degenaracy factor for the nth energy level. 
These are given as gv = 

47r/3fi3/2 X V/{2nf p|, ^7 = 2 x 3 and gn = n{n + l)/2 



Further, /x* is the chemical potential associated with the i-th quark. We may note that no 
correct thermodynamic limit can be considered when there is harmonic confinement as the 
system does not have translational invariance. One has to include center of mass motion 
corrections to be able to obtain correct thermodynamic limit. Such corrections, though 
important for hadron spectroscopy [9] and nucleon - nucleon scattering [10] involving few 
particles, are expected to be small (~ 1/N, N is the number of particles) and are not included. 

As we shall be considering hybrid stars, we shall consider the equation of state at zero 
temperature and for quark matter, along with u and d quarks, we include strange quarks 
and electrons which could be produced at high densities through weak interactions and be 
in equilibrium. The chemical potentials for each species of fermion may be written down in 
terms of the two independent chemical potentials fiB and fiE, the baryon and electric charge 
chemical potentials respectively as 

IJ'u = ^IJ'B + '^l^E (6a) 
= ^/is - ^f^E (66) 

f^s = ^/^B - ^f^E (6c) 
l^e = -IJ-E (6d) 

Thermodynamic pressure for quark matter is then given as 

Pqm ^ ^ Pi (7) 

i=uA,s,e 



where the electron pressure is given as 



127r2 



(8) 



These chemical potentials are fixed from the electrical charge neutrality condition for a given 
baryon density. The baryon number density and the electric charge densities are given as 



(9) 



and 



pE = -Pu - -Pd - -ps - Pe, 



where the quark densities at zero temperature are 



(10) 



Pi = J2^iP'-K)9v9i9n/V 



(11) 



and the electron density is 



37r2- 



(12) 



It is then straightforward to evaluate the sums in equation (5) for each quark fiavour to 
obtain the equation of state. 

In the neutron matter sector, the equation of state is still a matter of debate particularly 
at higher densities. Different parametrisations for the same are there for the description of 



the structure of neutron stars |jT3[. In contrast to sophisticated parametrisations designed 
for neutron star matter, we shall use throughout in this paper a simple parametrisation for 



the neutron matter equation of state by Sierk and Nix [|14 
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e7v(p) = P 



1) +Wo + Wsym + MN 



(13) 



9 M/Po 

In the above, Wq = — 16 MeV is the binding energy per nucleon at normal nuclear matter 
density, po=0.145 fm~^, Mjv=939 MeV, the rest mass of neutron and 14^5^^=32 MeV is the 
symmetry energy of neutron matter as estimated from liquid drop model calculations. While 



deriving the equation of state for the neutron matter, we take the value of compressibihty, 
K=550 MeV [0]. Below the nuclear matter densities the star matter is no longer determined 
by nucleon nucleon interactions only and we shall use for subnuclear density regions the 
equation of state given by Baym, Bethe and Pethick in the region 0.00267 fm~^ < 
p <1.6fm~^ and the equation of state given by Baym, Pethick and Sutherland for 



densities less than the above region. Different density regions of each equation of state were 
so taken that the pressure and energy density varied smoothly from one equation of state 
to the other PJl^ in the neutron matter sector. 



Once the energy density is known the pressure at zero temperature can be obtained from 
the thermodynamic relation |jl5| 



Pnm{.p) = - e^V, (14) 

and the chemical potential /ie is given as 

= (15) 

We shall now consider the scenario of phase transition from cold neutron matter to 
quark matter. As usual, the phase boundary of the coexistance region will be given by 
Gibbs criteria. The critical pressure and the critical chemical potentials are given by the 
condition 

PnmifJ'B) = PqmifJ'B)- (16) 

The left hand side of the above equation is given through equations (13) and (14) and the 
right hand side is the zero temperature limit of equations (5) and (7). 

We may note that in the calculations of pressure and chemical potentials in equations 
(5), (9) and (10) in the quark matter sector the zero point energy and in equations (13), 
(14) and (15) for the neutron matter sector the nucleon mass have been subtracted. This is 
because we wish to compare the pressure and chemical potential to consider phase transition 
for which the ground state energy is not relevant. However these will be included in the 
energy densities when we consider hybrid star as they contribute to the gravitational effects. 



In Fig. 1, we plot (P,/i_B) curves for quark matter and neutron matter which yield the 
critical parameters Per = 300MeV"//m^ and (/iB)cr=610 MeV. In Fig. 2, we have plotted 
the baryon number density in the quark matter sector as a function of the baryon chemical 
potential. The shape of the curve is indicative of the discrete energy levels characteristic of an 
oscillator potential. The baryon number density in the quark matter sector corresponding to 
the critical /i^ is calculated to be = 1.341/m~^. The same in the neutron matter sector 
at the critical is PB'"=0.784fm~^, about 5.4 times the nuclear matter density, which seems 
reasonable. At the critical pressure, the energy densities for the quark matter and neutron 
matter sectors, e^™ and e"^ are found to be 2542 MeV//m^ and 918 MeV//m'^ respectively. 
The discontinuity in the number density as well as energy density indicates a first order 
phase transition. We may note here that the calculation for the two phases are done in two 
different models-a phenomenological parametrisation in the neutron matter sector and the 
oscillator confinement model in the quark matter sector. The different treatments of the 
two phases essentially leads to a first order phase transition as above. 



For the description of neutron star, which is highly concentrated matter so that the 
metric of space-time geometry is curved, one has to apply Einstein's general theory of 
relativity. The space-time geometry of a spherical neutron star described by a metric which 



The equations which determine the star structure and the geometry are, in dimensionless 
forms. 



III. HYBRID STARS 



in Schwarzschild coordinates has the form 0,|TB|] 



ds^ = - 



e^^'Uf + [1 - 2M{r)/r]~^dr'^ + r^ldQ"^ + sm^Bc/^^] 



(17) 



dPirro) 



-G 



[e(fro) + -P(fro)] [M(rro) + 47raf^P(rro)] 
r2[l - 2GM{fro)/r] 



(18) 



df 




(19) 
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and the metric function, z/(r), relating the element of time at r = oo is given by 

di^jfrp) ^ ^^ [M(fro) + 47raf^P(fro)] 
df ~ f^[l-2GM{fro)/f] 

In equations (18) to (20) the following substitutions have been made. 

e = e/e„ P = P/e,, f = r/ro, M = M/M^, (21) 

where, with fi = 197.329 MeV fm and ro = 3 x 10^^ fm, we have 

a = e,rl/MQ, G = {G / h) / {r,/ M,,) (22) 

In the above, quantities with hats are dimentionless. G in equation (^2]) denotes the gravi- 
tational constant (G = 6.707934 x 10"^^ MeV'^). 

Further the surface condition for the metric function is given by 

v{R) = ln{l - 2GM/R) 

In order to construct a stellar model, one has to integrate equations (ff3) to (EUp from 



the star's center at r=0 with a given central energy density ec as input until the pressure 
P(r) at the surface vanishes. As stated, with central density greater than the critical energy 
density as dictated by the last section, we expect that we shall have quark matter, and not 
neutron matter at the center of the star. Hence we shall be using here the equation of state 
for quark matter using the zero temperature limit of equation (7) with -P(O) = -P(ec). We 
then integrate the TOV equations until the pressure and density decrease to their critical 
values, so that there is a first order phase transition from quark matter to neutron matter 
at radius r = rc for a given central density. For r > Tc, we shall have equation of state for 
neutron matter where pressure will change continuously (but the energy density will have 
a discontinuity at r = Tc). The TOV equations with neutron matter equation of state as 
given in equations (13) and (14) are continued. As described earlier we use the equation of 
state by Baym, Bethe and Pethick [O for subnuclear density region of neutron rich nuclei 



and of Baym, Pethick and Sutherland |]T8| for the lower density crystalline lattice bathed 



in relativistic electrons. The TOV equations are continued with the above equations of 
state till the pressure vanishes which defines the surface of the star. This will complete the 
calculations for stellar model for a hybrid "neutron" star, whose mass and radius can be 
calculated for different central densities. 

The energy density profile obtained from (|TB|) to (]2^) are plotted in Fig. 3 for central 
densities tc = 3250 MeV/fm'^ and ec = 4000 MeV/fm^. For core energy densities greater than 
the critical energy density (~ 2542 MeV/fm^) the core consists of quark matter. As we go 
away from the core towards the surface through TOV equations, when the critical pressure 
is reached, the density drops discontiniously indicating a first order phase transition. Thus 
e.g. for central density of 3250 MeV/fm^ such a star has a quark matter core of radius 1 km 
with a crust of neutron matter of about 10 kms, whereas for ec=4000 MeV//m^, the quark 
matter core radius is l.Skms with neutron matter crust of 9 kms. 

Fig. 4 depicts mass of a hybrid star as a function of its central energy density. Here we 
shall have two branches of solutions. Pure neutron matter stars at lower central densities 
Cc < e^*^ and hybrid stars at central energy densities ec > e^^. Taking into account the 
stability of such stars under density fiuctuations require dM/dec > PD[. We find stable 



hybrid stars with central energy densities Cc ~ 3 — A.GGeV/ fm^ and radii around 10 to 
11 kilometers, the critical mass being around 1.98 Mq, beyond which they may collapse to 
black holes. The instability region is shown as dashed curves in the above figure. 

Fig. 5 describes the mass radius relationship for hybrid star. This is important to 
estimate how rapidly the stars with a given equation of state can rotate without mass 
shedding at the equator. 



We may take the relativistic Keplerian angular velocity ^Ik given as pi 

0.72, 



M/Mq 



(23) 



lO^sec-i ■ \ (R/lOkm)^' 

as for neutron stars and estimate the same. Figure 6 shows the variation of Qk as a function 
of its mass for such a star. 

We wish to estimate the moment of inertia of these stars. The expression for the moment 
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of inertia of pulsars as well as the equillibrium for the rotating hybrid stars is given as 
. Stt 



^"■0 ,..4[e(^ro) + P(fro)] exp(-z/(rro)/2) 
arr ^ = 

1 - 2GM{fro)/r] 



(24) 



We further estimate the surface gravitational red shift Zg of photons as 



(25) 

[1 - 2G M/R] 

This is plotted in Figure 7 as a function of mass. In this context it may be worth while to 
mention here that the surface redshifts as determined from gamma ray bursters seem to lie 



in the range 0.2 to 0.5 |2^ where as the masses seem to lie in the range 1 to 1.85 Mq with 



an error of ±0.3 at both ends 26 



IV. CONCLUSIONS 

Let us summarise the findings of the present paper. 

Using a relativistic harmonic confinement model for quark matter and with a phe- 
nomenological parametrisation for neutron matter |T^ we saw that a first order phase transi- 
tion exists between the neutron phase and quark phase at about five to six times the nuclear 
matter density. 

We have used here the paramters of the harmonic confinement model fixed from hadronic 
properties. In this context or otherwise, we may note that the harmonic confinement model 
has its origin through Prasad - Sommerfeld dyon like configurations in QCD [§,0. When 
we include temperature and/or density effects we may expect that they will affect such 
configurations in a dynamical manner [0,^ and the parametrs of the confinement model 
may change. However to get an insight into the equation of state including the effect of 
harmonic confinement we have approximated Q to be independant of temperature and/or 
density. 

To study the possibility of a hybrid star, namely a star consisting of both quark matter 
and neutron matter, we applied the TOV equation to the appropriate equation of state with 
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a given central energy density ec- It turns out that a stable hybrid star with a quark core 
and a neutron matter crust can exist upto ec — 4.6 GeV/fm^ beyond which instability may 
result. For Cc from 3000 to 4600 MeV/fm^, the mass of the star varies between 1.795 to 1.98 
Mq and the radius of the star between 10 to 11 kms. Consistently, the bulk of the hybrid 
star is provided by the neutron matter, the quark matter providing core of about 1 to 2 
kms. The detailed properties of neutron matter e.g. whether it is soft or stif f do not seem 
to be reflected in the gross properties of hybrid stars. 

ACKNOWLEDGMENTS 

The authors are thankful to A.R. Prasanna, J.C. Parikh for discussions. AM and HM 
acknowledge many discussions with S.P. Misra, N. Barik and P.K. Panda. 



12 



REFERENCES 

[1] N. K. Glendenning and F. Weber, preprint LBL-31613, 1992; N.K. Glendenning, Phys. 
Rev. D46 (1992) 1274. 

[2] A. Rosenhauer et al Nucl. Phys. A540 (1992) 630. 

[3] H. Mishra, S.P. Misra, P.K. Panda and B.K. Parida, Int. J. Mod. Phys. E2 (1993) 547. 

[4] T. 0vergard and E. 0vergard, Class. Quantum Grav. 8 1191 L49. 

[5] T. 0vergard and E. 0vergard, Astro. Astrophys. 243, 412 (1991); ibid. Class. Quantum 
Grav. 8, L49, 1991; B.Datta, P. K. Sahu, J. D. Anand and A. Goyal, Phys. Lett. B283, 
313 (1992). 

[6] A. Mishra, H. Mishra, S.P. Misra and P.K. Panda preprint- IP/BBSR/94-1, PRL-TH- 
94/2, To appear in Z. Phys. C. 

[7] J. W. Negle, Nucl. Phys. A55 (1993) 47c; E. Manousakis and J. Polonyi, Phys. Rev. 
Lett. 58 (1987) 847. 

[8] S.B. Khadkikar, J.C. Parikh and P.C. Vinodkumar, Mod. Phys. Lett. A8 (1993) 749; P. 
C. Vinodkumar and S. B. Khadkikar, Phys. Lett. B329 (1994) 81. 

[9] S.B. Khadkikar and S.K. Gupta, Phys. Lett. B124 (1983) 523; S.B. Khadkikar, Pramana 
- J Phys. 24 (1985) 63; S.K. Gupta and S.B. Khadkikar, Phys. Rev. D36 (1987) 307. 

[10] S.B. Khadkikar and K.B. Vijaykumar Phys. Lett. B254 (1991) 320. 

[11] H. Huang, Statistical Mechnics (Wiley, 1963). 

[12] There is a difference of factor 47r/3 as compared to Ref. 7 in the expression of gv- When 
this factor is taken into account, the critical temperature of Ref. 7 reduces from 170 
MeV to 110 MeV. 

[13] W.D. Arnett and R. L. Bowers, Astrophys. J. Suppl. Series 33 (1977) 415. 



13 



[14] A.J. Sierk and J.R. Nix, Phys. Rev. C22 (1980) 1920. 

[15] See for example in " Quantum Theory of Many Particle System" A.L. Fetter and J.D. 
Walecka, McGraw Hill Book Company, 1971. 

[16] F. Weber and M.K. Weigel, Nucl. Phys. A493 (1989) 549. 

[17] G. Baym, H.A. Bethe and C. J. Pethick, Nucl. Phys. A175 (1971) 225. 

[18] G. Baym, C.J. Pethick and P. Sutherland, Ap. J. 170 (1971) 299. 

[19] S. L. Shapiro and S. A. Teukolsky, Black holes, white dwarfs and neutron stars (Wiley, 
New York, 1983). 

[20] S. Weinberg, Gravitation and cosmology (Wiley, New York, 1972). 

[21] J. L. Friedman, J.R. Ipser and L. Parker, Phys. Rev. Lett. 62 (1989) 3015. 

[22] K. Brecher, Astro. Phys. J. 215 (1977) L17. 

[23] J.B. Hartle, Phys. Reports 47 (1978) 201. 

[24] N.K. Glendenning, F. Weber and S.A. Moszkowski, Phys. Rev. C45 (1992) 844. 
[25] E.P. Liang Astrophys. J. 304 (1986) 682. 

[26] P.C. Joss and S.A. Rappaport, Annu. Rev. Astron. Astrophys. 22 (1984) 537. 

[27] A. Mishra, H. Mishra, S.P. Misra and S.N. Nayak, Z. Phys. C57 (1993) 233; A. Mishra, 
H. Mishra and S.P. Misra, Z. Phys. C59 (1993) 159. 



14 



Figure captions 

Fig.l: We plot pressure, P in MeV/£m^ as a function of the baryon chemical potential, hb 
in MeV for the quark matter (the sohd line) and neutron matter (dashed hne) . 

Fig. 2: We plot here baryon number density, in fm~^ as a function of baryon chemical 
potential //b in MeV for quark matter. 

Fig. 3: We plot here the energy profile curves inside the hybrid star for central densities 
3250 MeV/fm^ (the solid hne) and 4000 MeV/fm^ (dot dashed line). Discontinuity at the 
critical energy density is shown by the dashed line. 

Fig.4: Mass of the hybrid star as a function of central density, ec in MeV/fm^ is plotted 
here. The dashed portion indicates instability. 

Fig. 5: We plot here mass of the star in units of solar mass as a function of radius of the 
star in kilometers. 

Fig.6: We plot here the Keplerian angular velocity, in lO^sec"^ as a function of M/Mq. 
Fig. 7: We plot here the surface gravitational red shift, Zg as a function of M/Mq. 
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